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Experiments have been performed using microscopic beads 
to probe the smaU scale mechanics of actin solutions. We show 
that that there are a number of regimes possible as a function 
of the size of the probing particle. In certain cases we ar- 
gue that the quasi-static response resembles a smectic crystal 
rather than an isotropic solid, implying an anomalous scaling 
of the mechanical response of actin solutions as a function of 
the size of the probing particles. 
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The mechanics and rheology of actin filaments are a 
beautiful model system for the study of the dynamics 
and mechanics of semi-dilute polymers 0J^. They are 
characterized by length scales which are easily accessi- 
ble with optical techniques allowing the detailed study of 
phenomena such as tube dynamics. However, the macro- 
scopic rheology of these systems has been hard to master 
from the experimental point of view. Difficulties of purifi- 
cation and sample preparation lead to orders of magni- 
tude variations in such fundamental objects such as the 
value of the plateau modulus the standard mea- 

sure of the response of an entangled polymer solution to 
external perturbations. 

To get around these problems of macroscopic sam- 
ple preparation and also to probe the local viscoelastic 
behaviour of these materials a number of experimental 
groups have started using small, colloidal beads to study 
the local mechanics of these materials [pj-p^ . One cither 
pulls on the particles using super-paramagnetic beads in 
a magnetic field, or one simply observes the fluctuations 
of the particles undergoing Brownian motion. In this let- 
ter I shall try to attack the problem as to what exactly 
one measures in these experiments. In particular how 
large do these particles have to be in order to measure a 
macroscopic elastic modulus and when do we expect to 
be sensitive to the individual filament properties? 

In contrast with fiexible polymers solutions there are 
two principal length scales present in a semi-dilute solu- 
tion of actin: the mesh size and the persistence length. 
Naive application of scaling ideas thus becomes a highly 
ambiguous exercise because an arbitrarily large number 
of intermediated lengths can be created by considering 
with ^ the mesh size and Ip the persistence length. 
This ambiguity in lengths also translates into an ambi- 
guity in the plateau modulus which can be expressed 
as fcsT per characteristic volume. As an example of 
this difficulty we might quote two attempts to calculate 



the modulus in actin solutions with scaling approaches 
]r^ , p^ where two completely different results are found 
due in part due to this problem. Indeed this proliferation 
of lengths is already known for the tube geometry where 
one finds both a = +1/5 and a = -1/5, [^6|-@. We 
shall show in this article that a new intermediate scale 
with a = 3/5 becomes crucial in the understanding of the 
elasticity of actin solutions at length scales probed with 
micrometer sized beads. At these scales we show that the 
elastic response is highly anisotropic and resembles that 
of a smectic, with anomalous penetration of the response 
into the sample and unusual scaling of the response with 
the size of the probing particle. 

Note that in this letter I am interested in the low fre- 
quency mechanics and thus I exclude from the discussion 
high frequency fluctuation measurements (up to 20 KHz) 
which have been recently performed and am inter- 
ested in a quasi-static regime between lO'^Hz and lOHz. 
The reason for considering this range scale will become 
clear during the discussion. 

A coherent picture of the large scale mechanics of non- 
crosslinked actin solutions is now available. The actin 
system is usually polymerized pd| | in conditions such 
that the mean distance between filaments, ^ is between 
0.3/i and ^ can be linked with the concentration of 
monomers c by noting that ^ ^ 1/Vcd with d the size 
of actin monomers. A useful geometric quantity is the 
length of filament per unit volume p ^ 1/C^- The fila- 
ment is characterized by its persistence length Ip which 
is close to 15/x For a single weakly bent filaments 

the energy of a configuration is given by 



E = kBTlp/2 {d1v^{s)fds 



(1) 



where r_L(s) is the transverse fiuctuation of the filament 
about its equilibrium shape. 

In a manner which is familiar from flexible polymers 
the individual filaments are confined to a tube whose di- 
ameter scales as S,^^^/lp^^ and the filament is confined 
to the tube by collisions between the filament and its 
neighbors every 1^ ~ £,^^^1^^ le is is in some 



ways equivalent to the entanglement length in the Doi- 
Edwards tube model pl|-p^. The long time dynamics 
and mechanics are dominated by the reptation of fila- 
ments along their tubes ||5|, |23| , ^ . This process has a 
characteristic time, the reptation time, which defines the 
time scale beyond which the sample behaves like a vis- 
cous fluid (rather than an elastic solid) and can be as long 
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as several hours g]. Under macroscopic shear the longi- 
tudinal stresses in a filament relax relatively rapidly p^ ] 
leaving a residual contribution to the free energy which 
comes from the modification of the free energy of con- 
finement of the filament in its tube. A simple argument 
for this free energy is to count fc^T per collision of the 
tube with the filament. Thus the macroscopic modulus 
varies as G ~ pksT /le ^ c^ '*/Zp^^ as confirmed by an 
explicit calculation ||2^ , |2^ . 

This picture of filaments confined to a tube is only 
true on time scales that are long enough for the filament 
to dynamically sample fluctuations on the scale of le- 
This time, which is determined by the bending elastic- 
ity of the filaments varies as Te ~ rjle/lpkBT ~ lOHz 
. This is our reason for restricting our treatment to 
lower frequencies, at higher frequencies one is presum- 
ably sensitive to individual filament dynamics (coupled 
by hydrodynamics) rather than the collective, entangled, 
modes that interest us in this letter. For frequencies lower 
than the inverse reptation time (ie frequencies compara- 
ble to 10~^Hz) the sample behaves as a fluid and the 
bead moves freely as filaments slide out of the way of the 
particles. 

Before passing to the problem of the behavior of actin 
solutions we shall revise a Peierls like argument from 
which we can deduce the basic scaling behaviour of a 
normal elastic solid. We shall then adapt this argument 
to the case of semiflexible filaments: Consider a bead of 
radius R embedded in an elastic medium in d dimensions. 
If we pull on the particle with a force / we can make the 
following variational ansatz in order to find the minimum 
energy configuration. Let us assume that the material is 
disturbed over a distance / from the bead then the elastic 
energy, will scale in the following manner 

E,ar J {VafdV (2) 

where a is an amplitude of displacement, G an elastic 
constant and the integral is over the variational volume 
V ^ l"^. This scales as 

E,ar - G{a/lff (3) 

We see that in less than two dimension an arbitrarily 
small force is able to displace the bead large distances 
because E^ar can be made small by increasing the vari- 
ational parameter I. In three dimensions, however, the 
energy diverges with / and has a lower bound for small I 
due to the short wavelength cutoff coming from the finite 
size of the bead. Thus the minimum energy is found for 
I ~ R and we deduce that Eyar ~ Ga^R. At constant 
force the displacement scales inversely with the bead size, 

a « f/GR (4) 

A full calculation of the response of an isotropic viscoelas- 
tic material has recently been performed and confirms 
this simple scaling argument [p2[. 



We see that there is an anisotropy in the problem com- 
ing from the direction in which we apply the force /, and 
we should worry that the volume excited is not spherical 
as has been assumed in the argument. Let us perform 
a slightly more elaborate variational treatment where we 
assume that the volume V is characterized by an disk of 
dimensions I x I x D where the particle excites modes of 
wavelength I which penetrate D into the sample in the 
direction of /. In this case our estimate for E^ar is 

Eyar - {l^D)G{{a/lf + {a/Df ) (5) 

Where a/l and a/D are the estimates of the components 
of the strain tensor in the material. Taking Z) as a vari- 
ational parameter one sees that D ~ I and the problem 
reduces to that considered above. This is in fact a crude 
statement of the principle of St. Vernet that a force on 
a body with a wavelength I decays into the body over 
the same length scale, which is a elementary property of 
periodic harmonic functions in three dimensions. 

How must this argument be modified in the actin sys- 
tem? Experiments are performed with bead which vary 
in size from and 10/^. The smallest beads pass be- 
tween the filaments and diffuse almost freely ijllj; they 
will not concern us any further. Are we able to use con- 
tinuum elastic arguments (like that above) to deduce the 
experimental stress stain relationships? We now argue 
that in actin solutions there are now two contributions 
to the variational energy Eyar- For large beads the nor- 
mal continuum elasticity (summarized above) dominates, 
for smaller beads however a new, and novel elastic re- 
sponse is found: Consider a volume V distorted by a 
force on a particle of size R- Again we take this volume 
as anisotropic with dimensions I x I x D- In this vol- 
ume the filaments which traverse the volume bend with 
a wavelength I and there is a bending contribution to the 
total energy, coming from eq. (|^) which varies as 

El ^ {fD)ia^kBTlp/l^)p (6) 

The three multiplicative factors are respectively the vol- 
ume excited, the bending energy per unit length of fil- 
ament and the filament density within the volume, a 
is again the typical amplitude of the excitation in the 
volume. To this bending contribution one must add the 
equivalent of Eyar - When we impose the bending on the 
volume V there is also a variation in the geometry of the 
confining tubes. For instance in the direction of / the 
tubes are compressed by a factor comparable to a/D. 
Thus there is thus a contribution to the energy coming 
from the macroscopic bulk modulus of the form 

E2 - {l^D){{a/lf + {a/Df){pkBT/le) (7) 

where we have again respectively the volume, the square 
elastic stress and the macroscopic elastic modulus. We 
can now optimize Ei -f E2 by minimizing over D- How- 
ever, we first notice that there are two term linear in D 
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and that depending on the value of / one or the other will 
dominate. If we look at wavelengths 

; < ?c - yu; - (8) 

the contribution from Ei dominates over that from £'2- 
When I > Ic the second contribution dominates. 

We conclude that there is an important new length 
scale in the problem. When we look at excitations with 
wavelengths greater than Ic the two contributions in E2 
are going to dominate the elasticity and we are back to 
the case of normal continuum elastic theory. However in 
the short wavelength limit I < ^/Tj^ we find that the 
elastic energy is given by 

Eeff ~ a^p{fD){kBT/DX + (9) 

Minimizing the energy over D we find several surprising 
results. Firstly 

f ^ D^JTTp (10) 

the theorem of St. Vernet does not apply. Secondly 
substituting eq. ( [lO| ) in gives 

E,ff ^ pa^ksT^^ (11) 

which should be compared with the corresponding result 
for a normal solid eq. (||). From eq. (|l^) we see that 
the volume excited scales in an anisotropic fashion with 
the wavelength, quite unlike normal elastic solids. One is 
reminded of the penetration of excitation into a smectic 
liquid crystal with an effective energy for fluctuations of 
the form E ^ K J[{d^u)^ + (d^u)^ + P{d,u)^]dV where 
the z axis is defined by the direction of application of the 
force. 

The length scale I is absent in the energy (|ll|) which 
has not been minimized over the wavelength; this is anal- 
ogous to the case of a normal elastic material eq. (||) in 
two rather than three dimensions. It suggests that we 
are in the lower critical dimension for the problem (at 
least in a certain range of wavelengths) and that a fuller 
treatment will bring out logarithmic corrections to our 
picture. We also deduce, since eq. ( [Tl| ) is independent of 
I, that the amplitude of displacement of a particle should 
be independent of its size, in contrast to the dependence 
discussed above eq. (Q): We are no longer dominated by 
the short wavelength cutoff in the energy integral eq. (|^) 
despite being in three dimensions 

Finally we note that for coherency in this picture we 
require that the depth of penetration of the excitation 
into the sample D be greater than the distance be- 
tween the filaments ^, otherwise a continuum descrip- 
tion as used here must break down. This implies that 

I > k = J i\Jlclp = Sj^'^^iy^^ . On wavelengths shorter 



than this, one is presumably sensitive to filaments di- 
rectly in contact with the probing particle and do not 
feel the three dimensional nature of the sample. 

Substituting typical values for material constants, ^ ~ 
.5/i, Ip ~ 15/i we find that le ~ l/J.. The crossover length 
scale Ic ^ 4:IjL. The short wavelength cutoff li ^ 1.5/i. We 
thus expect the following series of crossovers as a func- 
tion of probing wavelength, (a) For £^ < I < li one probes 
the bending of individual filaments, (b) For li < I < Ic 
collective excitations of the solution become important 
with anomalous penetration of the excitation into the 
sample, (c) For I > Ic the elasticity becomes isotropic. 
These crossovers are too closely spaced to be experimen- 
tally studied in great detail, however we conclude that to 
measure a valid macroscopic response function particle 
sizes should be substantially greater than Ic = 4/i. 

Until now we have only considered the low frequency 
response of a sample, that is for times long enough that 
all longitudinal stresses have relaxed along the tube. It 
has been shown that one expects two plateau moduli 
as a function of frequency. The low frequency plateau 
used in the above discussion comes from variation in 
tube geometry under sample shear. The second much 
larger contribution which dominates at higher frequen- 
cies comes from coupling of the shear to the longitudinal 
density fluctuations of the filament in its tube. Can we 
see the crossover between the low frequency and high 
frequency behaviour with micro-bead techniques? This 
question is difficult to answer, the static approach used 
above is not adapted to answer this dynamic question 
however we can certainly expect that the frequency of 
crossover between the two regimes will vary with the bead 
size. 

The regime of the high plateau in macroscopic rhe- 
ology is delimited by the two times Tc ^ 0.1s and 
Teilp/le)"^ ~ lOs. This second time is the time needed 
for excitations to diffuse a distance Ip along the tube. It 
is important because macroscopic shear produces density 
fluctuations along the tube which are coherent over a dis- 
tance Ip. When we excite a sample with a wavelength I, 
which is smaller the Ip, we expect that the window of 
times for the observation of this high plateau is reduced 
to the interval between Tg and Tc{l/le)'^- For the smallest 
beads this high second plateau should almost completely 
disappear. Even with larger beads the elastic modulus 
should be substantially underestimated over certain fre- 
quency ranges. More detailed discussion of this regime 
seems to be difficult without a detailed dynamic theory 
of the coupling of the bend and longitudinal degrees of 
freedom. 

To conclude actin mechanics shows a quite rich se- 
ries of crossover in the response function G{q,Lu). We 
have simple arguments for the wavevector dependence of 
this function at frequencies between 10~^Hz and lOHz. 
Further work requires a full dynamic theory of the cou- 
pling between bending and density fluctuations. G{q, uj) 
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is more complicated than might be expected; recent ex- 
periments which interpret elastic stress propagation in 
terms of an isotropic elastic theory characterized by an 
elastic constant and Poisson ratio may be missing 
some interesting physics on length scales smaller than 
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